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■ We address a recent theoretical discrepancy concerning the Kondo effect in quantum dots with an 
' even number of electrons where spin-singlet and -triplet states are nearly degenerate. We show that 

the discrepancy arises from the fact that the Kondo scaling involves many parameters, which makes 

■ the results depend on concrete microscopic models. We illustrate this by the scaling calculations of 
the Kondo temperature, Tk, as a function of the energy difference between the singlet and triplet 
states A. rK(A) decreases with increasing A, showing a crossover from a power law with a universal 

' ' exponent to that with a nonuniversal exponent. The crossover depends on the initial parameters of 

^ I the model. 

. 73.23.Hk, 72.f5.Qm, 85.35.Gv 

■ 'T' ' Recent observations of the Kondo effect in semiconductor quantum dots have attracted much attcntion.0 In quantum 
^ . dots, the number of electrons N is fixed by the Coulomb blockade to integer values and can be tuned by gate voltages. 

a The usual Kondo effect takes place for an odd TV. Here, the spin 1/2 is formed in the dot and it is coupled to the 
Fermi sea of conduction electrons in external leads through tunneling barriers. „ „ 

Strikingly the Kondo effect has also been observed for an even TV in some special cases.B D In sp-called "vertical" 
quantum dots, the ground state can be changed from a spin triplet to singlet by the magnetic fieldB Sasaki et al. have 
found a pronounced Kondo effect near the degeneracy point between the spin states. □ In this experiment t^e Zeeman 
splitting of the triplet state can be neglected owing to a small g factor in semiconductor heternstructuresJj 
. In our previous paper, we have suggested an explanation of this phenomenon within a modelE The Kondo temper- 
O ' ature Tk has been evaluated as a function of the energy difference between the singlet and triplet states. A, using the 
"poor man's" scaling method. We have shown that (i) T'k(A) is maximal around A = 0, (ii) for positive A, IklA) 
decreases with increasing A obeying a power law, rK(A) oc l/A''', and (iii) for negative A, the Kondo effect is not 
7— I [ relevant when |A| Tk(0). The exponent 7 is not universal (0 < 7 < 1) depending on a ratio of the initial coupling 
^ ■ constants. Our results have indicated an enhancement of the Kondo effect by the comnetition between singlet and 
I triplet states, qualitatively in agreement with the experimental results by Sasaki et,a^.B After our work, Pustilnik 
. and Glazman have considered a different model for this "triplet-singlet Kondo effect" .eI Their results are qualitatively 

■ the same as ours but quantitatively different: They have derived a power law of TkIA) with a universal exponent, 

O ' 7 = 2 + Vs. 

' The purpose of this brief report is to elucidate the discrepancy between the results mentioned above. It should 
, be noted that a microscopic model sets the initial values of coupling constants, which are different between in Refs. 
and 1^, whereas a system of the scaling equations describes the evolution of the coupling constants. If it were a 
one-parameter scaling, both models would result in the identical results. However, this is not the present case of 
multi-parameter scaling. The scaling equations are 



I 

X3 



o 



J(l)2 + ( J2 + j2)/2 
J(2)2 + (J2 + j2)/2 



(1) 



dJ^^^/dlnD = 

dJi/dlnD = -2//(j(i) + j(2))Ji + lyJ'Ji 

. dJz/dlnZ? = -2;/(j(i) + j(2))J2-z.J'J2, 

X' [ dJ'/dlnD = 8z^(Jf-j|), 

H : 

5t 1 when the energy scale D is much larger than |A|. This system of the equations encompasses five coupling constants 
{J^^\ J^'^\ Ji, J2, J') and exhibits a complicated non-linear behavior. Our previous results correspond to a unstable 
fixed line of these equations, whereas the results of Ref. |8|to a stable fixed line. We show that in general the coupling 
constants develop along the unstable fixed line at first and transit toward the stable fixed line with decreasing D. The 
stable fixed line is much less relevant than in the situation of one-parameter scaling. Indeed our numerical calculations 
indicate that it may take up to six decades to reach the neighborhood of the stable fixed line (see Fig. 1(a)). We find 
that Tk(A) exhibits a crossover between power laws with nonuniversal and universal exponents, which are determined 
by the unstable and stable fixed lines, respectively. Besides, we show that our previous calculations along the unstable 
fixed line provide a lower limit of Ik (A) in an analytical form. 
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In a Coulomb blockade region with an even N, we consider two electrons at the background of a singlet state of the 
other TV — 2 electrons. We regard the background as the vacuum |0). The two electrons occupy two levels of different 
orbital symmetry, £i,£2, and make a spin-triplet state {S ~ 1) or -singlet state {S = 0): \SAI) = = c?||d2tl0)i 

|10) = (l/\/2)(di^4i +^u4t)|0)' |1 - 1) = 4i4il0>' |00> = (l/V2)(Ci4^4^ -C24t4i)|0>' w^ere d}^ creates 
an electron with spin cr in level i. The energy difference between the spin states is denoted by A = Es=o — Es=i. 
The coefficients in the singlet state, Ci and C2 (|Cip 4- |C'2p — 2), should be determined by thcj-ielectron-electron 
interaction and one-electron level spacing 6 — 62 — ei. We have set Ci = C2 in our previous work,Ll whereas C2 — 
in Ref. ||. Now we investigate general cases with respect to Ci and C2. The other singlet states with higher energies 
are disregarded. 

The level i (= 1,2) in the dot is connected to two external leads, L, R, by the tunnel couphiig, VL,i, Va^i. We 
assume two channels in the leads; the orbital symmetry is conserved in the tunneling processes^ We perform the 



unitary transformation for the conduction electrons in the two leads,0 



{VI 



(0 



VL,iC^fik^)/Vi with Vi = y/\VL,i\'^ + |VR,ip, where c]^'i^^ is the annihilation operator of an electron in 



ii) 



lead a, with momentum k, spin a, and orbital symmetry i (= 1, 2). The mode cj^^ is coupled to the dot with Vi. The 



mode c^^ is decoupled from the lead, which shall be disregarded. The Hamiltonian of the leads and that of tunneling 
processes are written as 
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/ . ^k(j ^kcr^ 
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H.c), 



(3) 



respectively. The density of states v in the leads remains constant in the energy band of [—D, D]. 

In the Coulomb blockade region, the addition and extraction energies, = E{Ni: 1) — E{N) =F are much larger 
than the level broadening Ti = 'kvV^ [i = 1,2) and temperature T, where E{N) is the energy with A'' electrons in 
the dot and fi is the Fermi energy in the leads. We also assume that E^ ^ 1^1 7 In the similar way to Ref. |^, we 
integrate out the states with ± 1 electrons and obtain the effective low-energy Hamiltonian, 
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dot 



H 
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The Hamiltonian of the dot, -ffdot, reads 



-ffdot — ^ Esfl^ifsM, 



(5) 



S,M 



SM 



using pseudofermion operators /^^,^ (/sm) which create (annihilate) the state \SM). It is required that fsAifs 
1. H^^^ describes the spin-flip processes among three components of the triplet state. 
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The exchange couplings are j(') ^ V^/{2Ec) {i = 1,2), with l/E^ = 1/E+ + 1/E-, which are accompanied by the 
scattering of a conduction electron of channel i. ffS^i^o represents the conversion between the triplet and singlet 
states, with the interchannel scattering of a conduction electron. 



kk' 



-(Ji/iVoo + J2/oVio)(4V/cf/ - 4!f 4'/) + (1^2) 



(7) 



where Ji — CiViV2/{2Ec). The last term of HcS represents the scattering processes without spin- flip in the dot. 
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kk'a i=l,2 



■J ^k'a'^ka JimJIM + J ' C^'^ C^<^/oo/00 
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We calculate Tk using the poor man's scaling method. Q With changing the energy scale (bandwidth of the con- 
duction electrons) from D to D — \dD\, we renormalize the exchange couplings not to change the low-energy physics, 
within the second-order perturbation with respect to H^^^ + jjS=i^o _|_ i/^^. This procedure yields the scaling 
equations in two limits. For D 3> |A|, H^ot can be safely disregarded in Hcs- The scaling equations are given by Eqs. 
(0) with J' = J'(i) - J'(2) - J"(i) + J"(2) . For D < A, the ground state of the dot is a spin triplet and the singlet 
state can be disregarded. Then J^^^ and J^^^ evolve independently, 

= -2i/jW^ (9) 



d\nD 



whereas the other coupling constants do not change. The Kondo temperature is determined as the energy scale at 
which the coupling constants become so large that the perturbation breaks down. 

In our previous model|3 Ji = and the scaling equations (|l|) are closed without J'. The ratio of J^^^J^^^ is 
also fixed in the scaling by Eqs. We obtain the Kondo temperature as a function of A, as follows, (i) When 
|A| <C T'k(O), the scaling equations (|l|) remain valid till the scaling ends. This yields 

Tk(0) = i?oCxp[-l/2;.(j(i) + J(2))], (10) 

where Dq is the initial bandwidth, which is given by ViJ+i?^ Jll (ii) When A ^^-Poi the scaling equations (^ work in 
the whole scaling region. Tk is identical to that of a localized spin with 5 = ljl3 

rK((») =I?oexphl/2i.j(i)], (11) 

when j(^) > j(2). (iii) In the intermediate region of Tk(0) <C A ^ Dq, we match the solutions of Eqs. (|l|) and (^) at 
Z) ~ A and obtain a power law 

rK(A) = rK(o)-(rK(o)/A)^ (12) 

The exponent is given by 7 = J^^^/j'^^. (iv) For A < 0, Tk drops to zero suddenly at A ~ — Tk(0). 

In general situations of Ci 7^ C2, we find that the line of j'^'^^ / j'^^'^ — const, and Ji = J2 is not stable in the 
scaling by Eqs. (|l]). The renormalization flow follows this line at the beginning, but finally goes to a fixed point of 

= j(2) = 00, J(2)/J(i) ^ 1, Ji/J(i) = ^^2(2 + ^5), J2/J^^^ = 0, and J'/J^^) = -2(1 + ^/5), as the energy scale 
D decreases to Tk- Around this fixed point, the stable fixed line discussed in Ref. || is relevant. By expanding the 
coupling constants around the fixed point to the first ordcLof l/ln(Z?/rK), Pustilnik and Glazman have obtained a 
power law of Ik (A) with a universal exponent, 7 = 2 + a/s.eI However, the range of energy scale is very limited where 
the physical properties are determined by the stable fixed line, due to the complexity of the multi-parameter scaling. 

To elucidate this multi-parameter scaling, we solve the scaling equations numerically when A > 0; Eqs. (jl|) for 
D > A and Eqs. (^ for < A. The Kondo temperature Tk is determined as the energy scale D at which 
i'(J^^-' -I- J^^^) = 1. First we examine the scaling by Eqs. (Q), assuming A = 0. In Fig. 1(a), we show the ratio of 
(J(i) - j(2))/(j(i) + j(2)) as a function of \ogD/TK{0), whem Tk(0) is the Kondo temperature at A = 0. We set 
j(2)^j(i) I- 5- a,nd c, 0.3 in the initial condition.^ In the case of Ci — C2 (solid lines), the ratio is fixed 

at a nonuniversal value determined by the given condition. In the cases of C2 = (broken lines) and C2/C1 = 0.6 
(dotted lines), the ratio gradually deviates from the initial value with decreasing T), and finally goes to zero. This 
limit corresponds to the fixed point of the renormalization, j(2)/j(i) = 1. It may need to reduce the energy scale by 
six decades to reach its vicinity. (The case a with J(i) = J(2) is an exception, where the fixed point is on the line of 
J(2)/J(i) =const.) Similarly Fig. 1(b) indicates (Ji - J2)/(J^^^ + J^^^). In the case of Ci = C2, the ratio is fixed at 

zero (Ji = J2). Otherwise, the ratio goes to a value of the fixed point, \/2 + VS/ V^, as D decreases. 

Figure 2 presents the Kondo temperature as a function of A on a log-log scale; a, 71^2)/ j(i) = 1 and 6, 0.3. Both 
Tk and A are normalized by Tk(0). In all the cases, the Kondo temperature increases with decreasing A. In the case 
of Ci = C2 (solid fines), Tk(A)/Tk(0) obeys a power law, Eq. ([l2|), with a nonuniversal exponent, 7 = j'2)/j(i). 
In the cases of C2 = (broken lines) and C2/C1 = 0.6 (dotted lines), Tk(A) shows a crossover from the power law 
with 7 = j(2)/j(i) at large A/Tk(0) to that with 7 = 2 -I- ^5 « 4.2 at smafl A/Tk(0). The crossover depends on the 
initial values of C2/C1 and j(2) The universal exponent is seen in quite limited situations. 

Finally we compare the Kondo temperature with various C2/C1 and j(2)/j(i) in the initial condition. We choose 
j^( J^^-' + J^^^) = 0.1 which corresponds to the experimental situation with sa 5K and level broadening F « 1.5K.EI 
In Fig. 3, we show Tk in units of T>o exp[-l/2i^( J^^) + J^^^)] (which is identical to Tk(0) with Ci = C2, Eq. (|l^), as a 
function of A. We find that the case with Ci = C2 (solid lines) provides a lower limit of Tk for each value oiJ^'^^/j'^^^ 
(a, 1; &, 0.5; c, 0.3), whereas the case with C2 = (broken fines) gives a upper limit. The inset to Fig. 3 shows Tk(0) 
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as a function of (= tan2 6'). For fixed C2/C1, Tk(0) is maximal at = 1 {0/n = 0.25) and decreases 

with decreasing J'^^^J*^'^-'. In a case of C2 = and j'^^-'/J*-^-' = 1, the development of the coupling constants by Eqs. 
(|l|) is restricted in a subspace with J2{D) = and J'-^\D) = J'-^^D), and in consequence they reach the fixed point 
fastest. This results in the highest Tk for a given value of J^^-' + J^^\ As the initial values are deviated from the 
condition, the renormalization flow to the fixed point needs more "time" in a larger space. (If D reaches A before 
Tk, the coupling constants follow Eqs. (||) at Tk < D < A, which is common to all the cases.) Since our previous 
model with Ci = C2 is the farthest from the condition of C2 = 0, it yields a lower limit of the Kondo temperature. 

In conclusion, we have examined a generalized model for the Kondo effect which involves spin-singlet and -triplet 
states in quantum dots, to elucidate the discrepancy between the previous studies.0l3 The Kondo temperature Tk 
is calculated as a function of the energy difference between the states, A, using the scaling method. The function 
of Tk(A) shows a crossover between power laws with a nonuniversal exponent (7 = J^^-'/J^^^) ancl with a universal 
exponent (7 = 2 + \/5). Although our previous model with Ji = J2 (Ci = C2 in the singlet state)Q is not complete 
to discuss the triplet-singlet Kondo effect, it yields a lower limit of Tk(A) in an analytical form, Eqs. ([Til), and 
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Figure captions 

Fig. 1: The scaling of the coupling constants by Eqs. (|l|), assuming A = 0. (a) (J^^' — J^^^)/{J^^'> + J^^^) and (b) 
(ii - J2)/{J'-^^ + J'-^^) as functions of log i:'/rK(0), where Tk{0) is the Kondo temperature at A = 0. is a, 

1; 6, 0.5; and c, 0.3. The cases with Ci = C2 are drawn by solid hues, C2 = by broken lines, and C2/C1 = 0.6 by 
dotted lines. Note that three lines are overlapped for case a in (a). Initially, we choose z/(j(i) + J(2)) =0.01. 

Fig. 2: The Kondo temperature, Tk/TkCO), as a function of A/rK(0), on a log-log scale (A > 0). rK(0) is the 
Kondo temperature at A = 0. J^^^/j'^^ is a, 1 and b, 0.3. The cases with Ci = C2 are drawn by solid lines, C2 = 
by broken lines, and C2/C1 — 0.6 by dotted lines. Initially, we choose + J(2)) = 0.01. 

Fig. 3: The Kondo temperature Tk as a function of A (A > 0). The units of Tk is Dq exp[-l/2i^(j(i) + J*^))] 
where Dq is the initial bandwidth {Dq = V E+E^ w 5K in the experiment of Ref. ||). j'^'/J^^^ is a, 1; b, 0.5; and c, 
0.3. The cases with Ci = C2 are drawn by solid lines, C2 = by broken lines, and C2/C1 = 0.6 by dotted lines. As 
the initial condition, we choose + j'^^) ~ 0.1. Inset: The Kondo temperature at A = 0, Tk(0), as a function 

of 61, where tan^ 6* = 6l/7r = 0.25 (0.20, 0.15) for ^ 1 (0.5^ 0.3). 
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Figure 1 (M. Eto and Yu. V. Nazarov) 
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Figure 2 (M. Eto and Yu. V. Nazarov) 
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